Normal-Superfluid Interface for Polarized Fermion Gases 
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Recent experiments on imbalanced fermion gases have proved the existence of a sharp interface 
between a superfluid and a normal phase. We show that, at the lowest experimental temperatures, 
a temperature difference between N and SF phase can appear as a consequence of the blocking of 
energy transfer across the interface. Such blocking is a consequence of the existence of a SF gap, 
which causes low-energy normal particles to be reflected from the N-SF interface. Our quantitative 
analysis is based on the Hartree-Fock-Bogoliubov-de Gennes formalism, which allows us to give 
analytical expressions for the thermodynamic properties and characterize the possible interface 
scattering regimes, including the case of unequal masses. Our central result is that the thermal 
conductivity is exponentially small at the lowest experimental temperatures. 

PACS numbers: 03.75.Hh, 68.03.Cd, 68.08.Bc 



I. INTRODUCTION 

By mixing fermion species of unequal particle number, 
recent experiments have revealed many aspects of the ex- 
citing physics of imbalanced fermion gases [E S S, 0, S 
B these experiments, a trap is loaded with an ul- 

tracold spin mixture consisting of two lithium hyperfine 
states. The existence of an easily-accessible Feshbach 
resonance allows the tuning of the interspecies interac- 
tions with great accuracy so that the entire BEC-BCS 
crossover with number imbalance may be explored Q. 
The BCS regime involves a small and negative s-wave 
scattering length, with the fermions forming weakly 
bound Cooper pairs. On the other hand, in the BEC 
regime, where the scattering length is positive, the par- 
ticles form ti ghtly b ound molecular pairs forming a con- 
densate [1, [9I llCll.Tllj| . In between lies the unitary regime, 
characterized by a diverging scattering length (signify- 
ing the appearance of a bound state) , where particles are 
strongly interacting. The divergence of the scattering 
length implies that the system is universal, in the sense 
that only a single length scale characterizes it. 

Much work has been done in an attempt to understand 
the ground state properties of the imbalanced clouds; 
however, despite theoretical advances [H, [H, [3, [H, [H, 
[13, [3, [ll, 0, mi in M, some problems remain unre- 
solved. The MIT data El,i,i,[3 are in good agreement 
with Monte Carlo (MC) simulations [l^l and renormal- 
ization group (RG) predictions [l5|, which point to the 
existence of a three-shell structure consisting of a super- 
fluid (SF) core phase, a partially polarized and a fully- 
polarized normal (N) phase. From these experiments, the 
homogeneous phase diagram has been determined and is 
in agreement with both the MC and RG results. On 
the other hand, the results of the Rice group's experi- 
ments [H [3| are are less well understood. For example, 
they find that superfluidity does not break down at high 
polarizations, which is incompatible with the results of 
the MIT experiments and thus also in contradiction to 
theoretical predictions of a critical polarization, above 



which superfluidity breaks down. Since both the MIT 
phase diagram and the MC and RG predictions were 
obtained under the assumption that the local density 
approximation (LDA) holds, and since the Rice experi- 
ments involve smaller traps, it is likely that the answer to 
the riddle lies in the breakdown of LDA (in other words, 
finite-size effects); indeed, taking as an ansatz a two-shell 
structure, the Rice data can be recovered by using a N-SF 
interface tension as a fitting parameter [lO] . 

For our purposes, it is sufficient to say that the MIT 
experiments are well-described by theory, and furnish di- 
rect evidence that a N-SF interface exists, separating an 
unpolarized SF from a partially polarized N. 

In this article we argue that the N-SF interface may 
block the thermal equilibration process between the N 
and SF phases; one consequence of this will be to cause 
a temperature imbalance to appear across the interface. 
Incorporating such temperature difference in the existing 
models may provide the key to a better understanding 
of the experiments. We present a detailed study of the 
scattering off the N-SF interface within the Hartree-Fock- 
Bogoliubov-de Gennes (HFBdG) formalism, followed by 
a calculation of the heat conductivity across the interface. 
We find that this conductivity drops exponentially fast 
with increasing A/k^T with A the SF gap and T the 
temperature. Our model predicts the conductivity to be 
vanishingly small for the experiments at unitarity. A 
summary of our results was published earlier in Ref. [H. 

The article is structured as follows. We start off in 
section |lT] by discussing the bulk physics of imbalanced 
fermion gases. This is done using a mean-field Hamilto- 
nian which incorporates both the gap formation and the 
Hartree-Fock terms. We present analytical results for the 
equations of state which enable us to write down the con- 
dition for coexistence between the N and the SF phases. 
Then, in section IIIH we discuss in detail how particles 
scatter on a N-SF interface according to the Hartree- 
Fock-Bogoliubov-de Gennes formalism. We find a rich 
variety of interfacial scattering processes, depending on 
the energy and momentum of the incident particle. The 
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implications of this scattering are outlined in section ITV] 
where we obtain the thermal conductivity exactly and 
with a good analytical approximation. Finally in sec- 
tion |V] we calculate that for the experiments at lowest 
temperatures, the typical timescale associated with the 
equilibration of a temperature difference across the inter- 
face, is of order seconds and therefore comparable with 
the trap lifetime. 



II. BULK PHYSICS AND COEXISTENCE 
A. Bulk Phases 

Fermion mixtures in the BEC-BCS crossover are the 
subject of a very active research field. Although some 
theoretical models give accurate results for the BEC 
and/or BCS regimes, no model yet allows an analyti- 
cal, quantitative description across the entire crossover. 
In the following, we use the so-called one-channel model 
at zero temperature; this is a straightforward general- 
ization of the BCS model to arbitrarily strong coupling 
and is widely used in the literature, as it is known to em- 
body most qualitative features of the BEC-BCS crossover 
for the currently experimentally relevant Feshbach reso- 
nances. Interactions between fermions and the bosonic 
bound states may be explicitly incorporated within a 
two-channel model 9']. 

We start off by deriving the Hartree-Fock-Bogoliubov- 
de Gennes (HFBdG) equations. We then focus on 
the physics of an interface separating two semi-infinite 
phases: a SF and a polarized N. We write down the equa- 
tions of state in these phases, and use them to find the 
condition for two-phase coexistence. In this section, we 
extend the work reported in Ref. [2^ by including the HF 
term into the BdG equations and treating the case of 
unequal species masses. 

The system under consideration consists of two 
fermionic species t and | with masses rrii and chemi- 
cal potentials for i —1, i and which we assume to be 
trapped by the external potential Vi{r). The intraspecies 
interactions are via p-wave scattering which can be safely 
ignored at low temperatures. On the other hand, the in- 
teractions between the two spin species are taken to be 
contact interactions characterized by the coupling con- 



stant G = AttH a/m^ < with = 2(m^ 
and a the s-wave scattering length. Within a general- 
ized HF approximation, one may work with the following 
Hamiltonian 123, Hi : 



H = dr 



$|(r)WT(r)$T(r) + ^{r)Hi{r)^ i{r) (1) 
+A(r)$|(r)$j(r) + A*(r)$i(r)$t(r) 



where = -h'^V'^ /{2m,) - ^, + ^^(r) + J7,(r), the 
gap A = -G'($f (r)$|(r)) and the HF terms U-^ir) = 
G($|(r)$|(r)) and Ui{r) = G($|(r)$|(r)). As usual, 



one diagonalizes the Hamiltonian ((T]) by a Bogoliubov 
transformation i.e. by writing the operators '^i in terms 
of fermionic quasiparticle operators c: 



-iEit/n 



where Ui^i and vij are the quasiparticle wave functions, 
{ceA,ce'j} = and {cl,,Ci>j} = 6e^e'S,j. Then, if un- 
der this transformation the Hamiltonian ([T]) is diagonal- 
ized, the relations [H, cg^] = —E^c^^i and [H, cj ■] — 
Enc\ j must hold such that one straightforwardly finds 
the four Hartree-Fock-Bogoliubov-de Gennes (HFBdG) , 
or Blonder- Tinkham-Klapwijk, equations [25l . |27| : 



A*u£,-|- - Tiive^i = El 
A*ue^l - ^T^^,T = E!ivi^^ 



(2a) 
(2b) 
(2c) 
(2d) 



The wave functions j and vi^i which appear in the first 
two equations are not coupled to the wave functions u^j 
and of the last two equations so that, apart from 
normalization, they can be treated separately. 

In the rest of this section, we will be concerned with the 
study of a homogeneous system of two fermion species at 
fixed chemical potentials /i| and /i^ and with no external 
potential V-\ — Vi = 0. Furthermore we treat both the 
case of a N and a SF phase. 

At zero temperature and for negative scattering length, 
there always exists a (meta-)stable state with A = A* 7^ 
and for which both species densities are equal; we call 
this state the BGS state or SF phase. For this state, we 
define now the effective chemical potential imbalance h^ 
and the average effective chemical potential fi^ as follows: 



(3a) 
(3b) 



Without loss of generality, we further assume that > 
0. Note that in the SF phase, both spin species have 
equal densities which implies equal HF potentials = 
Ui = Ui for the spin species. One can then find the 
ground state energy of the homogeneous SF system as a 
function of fj,^, /Xj and a. Solving the HFBdG Eqs. 
and minimizing the ground state energy with respect to 
the variables A and i7g, the gap and number equations 
at zero temperature are [59| : 



1 = 




1 



v/(£,-Ms)^ + A2 s 



(4a) 
(4b) 
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where = k /(2m^). Note that the above equations 
are coupled by the term which appears in the defi- 
nition of . The ground state grand potential per unit 
volume in terms of A and /ig is: 



d^k 
(2^ 



AV2+(e 



(5) 



We are mostly interested in the BCS side of the BEC- 
BCS crossover where the effective chemical potentials 
fi- — \J. and fig can be taken as positive. Performing the 
integrals over the three-dimensional wave vectors leads to 
the analytical expressions for the gap, the HF potential 
and the grand potential per unit volume [H, [33, [H, [3^ : 



(TO^)3/2e5/2r,2 
"80^fi3|P,,°(7?)|5 



(6a) 
(6b) 



(6c) 



with P the Legendre function, r/ = — 1/ -^/l + (A/ii^)^ 
and we introduced the energy scale which is set by the 
scattering length = ti^ / (m^a^). Expression ([Sc]) al- 
lows to write the grand potential in terms of A/e^; in- 
deed, for each value thereof, one can find r/ by solving the 
relation: 



(7) 



such that the grand potential ([6c| is fully determined by 
A/e^ only. 

For the normal state, one finds the following expres- 
sions for the grand potential per unit volume of the 
ground state and for the HF potentials: 



El 

Ul 



2V2 



m^/'(^l -f/|)5/2 (8a) 



4%/2 /m^\ 
37r J 

4^/2 



37r 



3/2 



3/2 



Mi - Ul 



/i| - Ul 



3/2 



3/2 



(8b) 
(8c) 



At fixed chemical potentials, the attractive interactions 
induce an increased normal state density and a decreased 
ground-state grand potential. Eqs. (j8bp and ((8c|) are cou- 
pled in a non-trivial manner and can be solved numeri- 
cally to find [/| and C/j^. 



B. Bulk Two-Phase Coexistence 



Before writing down the N-SF coexistence condition, 
we explain how a variation of the chemical potential im- 
balance between the spin species may change the ground 
state. A small imbalance (or polarization) gives rise to 
a SF state while a large polarization leads to a N state. 
Intuitively this is expected since it is known from BCS 
theory that pair formation is energetically favorable for 
particles with opposite momenta. By shifting the Fermi 
energies, an energy cost must be paid which is associated 
with the matching of the Fermi levels so as to enable pair 
formation; when this energy cost is too high, the N state 
prevails. 

In a trapped imbalanced fermion system, both the 
N and the SF phase may be present, bounded by an 
interface. Using a sufficiently large number of parti- 
cles, a local density approximation is justified there and 
amounts to using an effective local chemical potential 
^i(r) = fx^ — Vi[Y) with M?_the chemical potential at the 

trap center and i =T, i [i,!!!!!^,!!!,!!!,!!^,!!!!- Since in 

experiments both spin species feel the same trapping po- 
tential, the imbalance m (r) — /ij^ (r) is constant through- 
out the trap. On the other hand, the average chemical 
potential [lJi]{v) -I- /i|(r)]/2 varies: it is large at the cen- 
ter of the trap and decreases towards the trap boundary. 
The ground state energy of the SF increases with increas- 
ing value of the gap; as seen from Eqs. (pa|) . the gap value 
is proportional to the average chemical potential /^g and 
independent of . Therefore, if present, the SF phase 
will be at the trap center where /z^ is maximal. Moving 
radially outwards, two-phase coexistence of the N and 
SF phase occurs at the position where there is a balance 
between the energy gained by creating, on the one hand, 
a gap in the SF, and, on the other hand, a density differ- 
ence in the N phase. 

To establish now the condition for coexistence between 
the N and the SF phase, recall that we have a system con- 
taining two fermion species i =t, i with masses rn; and 
chemical potentials /i^ for which interactions are charac- 
terized by a negative scattering length a. We aim to find 
the coexistence condition as a function of the parame- 
ters a, ^1 and The SF state is fully described by 
the parameters A and , which can be extracted from 
Eqs. ([6a|) . ()6b|) and The normal state, on the other 
hand, is characterized by and Ui which are solved for 
using Eqs. (l8b| and ([5c)l . 

The coexistence condition obtained by equating the 
grand potentials of Eqs. (f6c|) and (l8a|) is: 
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V2h^ 



(10) 



The Clogston limit is exact in the BCS hmit, when <C 

Lastly, it is common to write the scattering length 
a in units of l/fc^, the latter defined by the relations 

find that in the SF 



k„a 



3n {P,,M/P./M-V) 
4 PIM 



1/3 



m 



3/2 / jj \ 5/2 



X 3/2 , 



Ui 



5/2 



(9) 



By condition ^ one can find the surface of coexistence 
in the space of the parameters , /Lt|, a). 

For one specific case it is easy to relate the gap value to 
the chemical potentials at coexistence; this occurs in case 
of equal masses = m| = TO| without incorporation 
of the HF term. The energy scale can then be elim- 
inated and by a numerical analysis, one finds to a good 
approximation the Clogston relation, originally derived 
for superconductors in a magnetic field [42J : 



k^T^ = h^kl/{2m^) = ?i2[37r2n,]2/3/(2m+); we 




FIG. 1: The N-SF interface (thick vertical line) with the spec- 
tra of the species t and I on the N side and the gapped a, /3 
spectra in the SF. The long dashed line cuts the spectra at 
particle-like (filled dots) and hole-like (empty dots) quasipar- 
ticle states, all having the same energy. The minimal values 
of the a and the (3 spectrum are 2(\/mA — e° ^)/(m -|- 1) 
where ej^ ^ = [±/ig (1-1- m) + (1 — m)]/2. An incident parti- 
cle of species t (curly arrow) with momentum and energy 
E can have up to four scattering channels: the Andreev re- 
flected fc^ hole, the specularly reflected —k^ particle and the 
transmitted hole-like —k'^ and particle-like k'^ states. 



III. INTERFACE SCATTERING 



A. Bulk Particles 

In this section, we explain in full detail how particles 
scatter from a N-SF interface for imbalanced fermion 
gases. By taking a simple model for the interface, we 
describe the different scattering regimes and we give the 
transmission coefficients for scattering across the inter- 
face; the latter will be used for calculating transport 
quantities in the next section. 

To model the N-SF interface, we start with a geometry 
wherein the z = plane separates the N phase at z < 
from the SF phase at z > or: 



A(r) = e(z)A, 



(11) 



with Q the Heaviside function. This approximation 
was also used by Andreev [i^ for studying the normal- 
superconductor interface. 

A particle of energy E incident on the interface may 
scatter to states which have the same energy and these 
are deduced from the spectra in the SF and in the N 
phase. In the N phase, we have for species i =^,i: 



E, 



21.2 



2m, 



We depict with full lines this spectrum in the left panel 
of Fig. [T] where one recognizes two approximately linear 
energy branches "f and | near the (unequal) Fermi sur- 
faces. In the SF on the other hand, one finds using the 
HFBdG equations, an a and a (3 spectrum: 



E 



fi2k2 

2m 



[2^ 



A2, (13) 



where m^ = 2(m^ "'^ ± m, ^) ^. These gapped spec- 



'i _ 
tra are depicted in the right panel of Fig. [1] The a 

branch is composed of particles of phase t and holes 

of phase | and vice versa for the /? spectrum. One 

can find that the minimal values attained by the a and 

P spectra are 2(-\/mA — ^)/(m -I- 1) where ^ = 

[±hg{l+m)+^^(l — m)]/2 with m = m|/m|, whence the 

lowest spectrum touches the E = axis or becomes un- 

gapped under the condition that VmA < e° . We assume 

a gapped spectrum. Since the gap function is constant 

in both the N and the SF phase, we can decompose the 

(12) quasiparticle wave functions Uf and v^, associated with 

the a-branch, into their Fourier components: 
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T - 



k,± 

k,± 



k, ± 



(14) 



Bps ' ^ 

k,± 



Here A and B are complex scattering amplitudes and 
the sub- and superscripts n, s, p and /i denote normal, 
superfluid, particle-like and hole-like, respectively. We 
split the vectors k into its components parallel to the 
wall k|| — {kx,ky,0) and its ^-component fc^, k^, k^ and 
k*^ perpendicular to the interface. At fixed energy, there 
exist eight quasiparticle states, four of which in the N 
and four in the SF phase. These states are represented 
by dots in Fig. [TJ filled dots are particle-like states and 
empty dots are hole-like states (the latter have oppo- 
site group and phase velocity). By use of the spectra of 
Eqs. (fT^ and one relates the wavenumbers /c^, k*^ 
and k^ to the wavenumber k^ of a particle in the N phase: 

K = (-^« + (MT -U^+ EJ/H^^ , (15a) 
fc, - (fc2 + 2mT(C/t-mC/x (15b) 

t2\l/2 



+Us{m-l)-2eJ/hy 
(/c2 + 2mT(C/T-C/,-x^)/fi')'^', (15c) 



where we defined: 

_ {E^+hs){l + rh) ^s(l-m) 



and fh = mi /mi. 



(16a) 
(16b) 



B. Scattering Diagram 

Let us now use these equations to try to understand 
the physics of the scattering processes at the interface. 
Consider for example a particle incident on the interface 
from the N side with wavenumber k^ (along the z-axis) 
and energy ; this is represented in Fig. [T] by a curly 
arrow. Whatever processes occur at the interface, energy 
must be conserved; but out of the seven remaining states 
in the N phase with the same energy, only four have a 
group velocity directed away from the interface, this im- 
plies that B^ "^ = A^'''_ = B^'^ — 0. The particle can now 
undergo either specular reflection to the state marked in 
Fig.[l]by — fcp, Andreev reflection to the state fc,^ or it can 
penetrate the SF as a particle- like (to k^ ) or as a hole-like 
quasiparticle (to — fc'*). Note that Andreev reflection is 
also of importance in other ultracold systems such as in 
fermionic traps without population imbalance j43j and in 
bosonic traps near a N-SF interface in the absence [3| 
or presence [i^ of a lattice. 
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The scattering amplitudes depend on the energy E^ of 
the incident particle, represented by the long dashed line 
in Fig. [TJ Consider for instance decreasing the energy of 
the incident particle to a value smaller than the lowest 
value of the a-branch; since no SF states are available, 
transmission to the SF is then not possible. Mathemati- 
cally this is marked by the wavenumber of one state in the 
a-branch becoming imaginary, indicating the presence of 
an evanescent wave in the SF. Furthermore, we can also 
raise the long dashed Hue of Fig. [H that is, increase the 
energy i?^ , to a value above the local maximum at fc^ = 
of the spectrum of species J, in the N phase; this will in- 
hibit Andreev reflection. 

The crossover of either of the k. , A:'* and k^ and the 

h ' 

wave vector k|| from real to imaginary (or vice versa) 
signifies a change in the scattering mechanism. Defining 
Xp = fi^fc^/2m|, such transitions are encountered when 
(see Eq. ps]) ): 



(17) 



The curves 1-4 separate the possible scattering regimes 
and are shown in Figs. [5] and [3] where the latter depicts a 
general (gerenic) diagram and the former assumes equal 
species masses and zero HF potentials. 

To better understand the various scattering regimes, 
we switch to working with the energy of Eq. (jl6ap 
instead of E^ . Recall that takes values in the interval 
[£° , oo [ where £^ — [/i^ (I-I-to)-!-//^ (l—m)]/2; clearly = 
E^ in case = and TO| = m^. If the energy of the 
particle incident on the interface is small enough or < 
VmA, there is insufficient energy to excite quasiparticles 
inside the SF; the particle is then fully reflected as a 
superposition of a particle and a hole. This regime is 
labelled I in Figs. [2] and [31 The situation where E^ < 

— /i| + [/| i.e. regime II is physically forbidden since the 
wave vector of the incident particle (i.e. k^ ) is imaginary. 

Suppose now that a particle is incoming with energy 
above the threshold value VmA. Although a particle 
in regime III appears to have sufficient energy for trans- 
mission, in undergoes complete reflection. Particles in 
regime IV may be transmitted and reflected but all tar- 
geted states are particle-like. In both regimes IV and V, 
there are no reflected holes; that is, Andreev reflection 
does not occur. However, in regime V, both particle-like 
and hole- like excitations are present in the SF. Finally, 
in regime VI, both particle-like and hole-like excitations, 
or both Andreev and normal reflection, are allowed. 

The blocking of scattering regimes can be understood 
as being caused by a too large angle of incidence at the 
interface, similar to total internal reflection. While the 
critical angle is roughly constant for optical total internal 
reflection, here the critical angles of incidence depend on 
the energy of the incoming particle. Generally, a particle 
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FIG. 2: The various scattering regimes for a particle of 
species t which is incident on the interface from the N 
side, as a function of its "energy" (see Eq. (I16a|l ') and 
Xp ~ h^k^/2mf (where is its momentum along the z- 
axis). We assumed here no species mass asymmetry m| = 
and no HF potentials (7m, s = 0; we depict a generic dia- 
gram for the more general case (m| 7^ m| and I7|,x,s 7^ 0) 
in Fig. [3l The heavily-shaded regimes are energetically for- 
bidden while complete reflection occurs in the lightly-shaded 
regimes. Particles in regime VI may scatter to all states in- 
dicated in Fig. [1] by curly arrows. Above line 4 (regimes IV 
and V), Andreev reflection does not occur, and above curve 
3 (regime IV), hole-like excitations are also impossible. The 
functions of the curves 1-4 are given in Eq. (|17p . 



of species t has an angle of incidence 9 which satisfies 
tan6' = |k l/fc , or in terms of the energy: 



tan(0)=('^^I^^i±^-l 



1/2 



(18) 



The critical angles for transmission, hole-like transmis- 
sion and Andreev reflection can then be obtained as a 
function of the energy by substitution of x from 
Eq. ini) into Eq. (HH). ° 

To summarize, the lightly-shaded regimes in Figs. [2] 
and [3] describe normal particles of species t which un- 
dergo complete reflection, while the rightmost, heavily- 
shaded regime is unphysical. Only incident particles in 
the unshaded regimes may excite quasiparticles inside the 
SF. If we consider now holes of spin species | which are 
incident on the interface from the N phase, they can only 
be transmitted in a regime analogous to VI in the dia- 
gram. Furthermore, one can flnd now also the scattering 
diagram, applied for particles of spin species | and holes 
of spin species t incident on the interface and scattering 
to the /3 branch. The scattering diagram is then quali- 
tatively similar to that of Fig. [2] and the curves 1-4 are 
given by Eqs. (jl7p after performing the transformations 
/13 —hs, rh 1/fh, a — > /3 and t^i- 



FIG. 3: The same applies as in Fig. [5] except here m| 7^ 
and f7|,|,s 7^ 0. Regime I is bound by threshold line — 

\/mA and the lower limit = e° = [hsi^ + ^) + Ms(l ~ 
m)]/2. Compared to Fig. [2] it can be seen that curves 1, 2, 3 
and 4 are shifted. 



C. Discussion of the Scattering Diagram 

At this stage, we can discuss the influence of the inter- 
action parameter k^a, the HF potentials and the species 
mass asymmetry on the diagram of Fig. [2] as they may 
strongly affect the transport properties which will be cal- 
culated in the following section. 

The interaction regime — First of all, the x "^i^" 
grams of the deep BCS interaction regime and the unitary 
interaction regime have a different topology. Assume for 
now that = and that ~ U]- = Ui = 0. In 
the deep BCS regime, the relation 2A ^ A + /i^ implies 
that line 2 in Fig.[2]is shifted horizontally to much higher 
values of x so that regime VI becomes by far the most 
important scattering regime. Therefore quasiparticle re- 
flections occur almost only via the Andreev mechanism. 
It is known that this involves transport of particles, but 
no transport of energy across the interface. In contrast, 
at unitarity where 2 A > A -I- /ig (since A ^ 1.16/j,j, as 



derived from Eq. (|6ap in case m-f = and t/|,x,s ~ 
phase VI does not even exist. This means that quasi- 
particles with energy above the threshold A cannot 
undergo Andreev reflection but only normal (specular) 
reflection, by which neither particles nor energy are car- 
ried across the interface. 

The HF potentials — Assume again that ~ rnj^; 
the HF potentials C/^, [/-f and are seen to be present 
in the expressions (|17p which describe the curves 1-4 of 
Fig. [21 However, it follows from these equations that the 
HF potentials induce only horizontal shifts of the curves 
1-4 but do not affect the upper and the lower bounds of 
regime I. Whereas line 2 is shifted to the right, the curves 
1, 3 and 4 are always shifted to the left. This can be 
deduced from the fact that Ug < < and < U-\ < 
since the interactions are attractive and the density of 
the SF exceeds the density of the N phase. 

The species mass asymmetry — In the presence of a 
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species mass asymmetry ^ m|, regime I may become 
either larger or smaller. This is due to a shift of its lower 
bounds" = [hg{l+m)+fig{l—rh)]/2 (obtained by taking 
= in Eq. (|16ap ) and a shift of the upper bound 
— VmA. It is seen that both shifts have the same 
sign but can either increase or decrease the surface area 
of regime I. Also, the species mass asymmetry induces a 
shift of line 4 as is seen from Eq. pT|) . 



A and B of Eqs. ((T4)) to the transmission coefficients. The 
amplitudes are determined by matching the quasiparticle 
wave functions of Eq. (|14p and their derivatives at the 
interface z — [46]. In its most general form, this leads 
to the equations: 



D. Transmission Coefficients 



In order to study the transport properties in the next 
section, we first seek to relate the scattering amplitudes 



I k,+ 

k,+ 



k, — 



\ ( 



k,+ k.+ 

q Bp-" + c B'^-" + (; 

^P k.+ k,+ 

V kpq (Bp'^ - Bp " + /fcH, (B^^.^ - ) / 



Ap-' + A'' 
k,- k 

Bp-^ 



-A" 



(19) 



where ±/A. As argued before, de- 

pending on the various scattering regimes, one must set 
some scattering amplitudes to zero. One finds: 

for regime VI: 
for regime V: 
for regime IV & V 
for holes: 

We give the solutions for the scattering amplitudes in 
expressions ((27|) -(p0 | in the Appendix. 

From the HFBdG Eqs. Q ft readily follows that the 
quasiparticle density p^{r) = |u|(r)p + |t;|(r)p and the 
quasiparticle currents 



k.- 


= Ap-^ 

k.- 


= 

k.+ 


= 0, 


k.+ 


= Ap-^ 

k.- 


k,+ 


= 0, 


k,+ 


= Ap-^ 

k,- 


k,- 


= 0, 


k,+ 


k,- 


k,+ 


= 0. 



2mi 



satisfy the "continuity equation" dp^/dt + V • =0. 
Note, however, that they do not express conservation of 
particles or mass. The same expressions are valid for 
the /3 branch after the transformation t^i- We define 
now the transmission coefficient Sp (e^ ,Xp) of an incident 
particle of energy and momentum k as the ratio of 
the transmitted to the incoming current along the z-axis. 
As argued before, the transmission coefficients vanish in 
regimes I, II and III of Fig. [21 for regimes IV, V and VI 
(and for holes) we give their analytical forms in expres- 
sion (|3ip in the Appendix. For energies slightly above 
the transmission threshold « VmA (for regimes V 
and VI) one can analytically find the following behavior 



[viVvl - vlVvi] 



This is similar to the case of particles scattering from a 
HE potential of height VmA ^ . Moreover, in the BCS 
regime wherein, by the Andreev approximation, one can 
assume kp = k*^ = k = k^, the transmission coefficients 

p h ' 

reduce to: 



Apart from the species mass asymmetry m, this is the 
expression which was obtained by Andreev in Ref. (48l |. 



IV. HEAT TRANSPORT THROUGH THE 
INTERFACE 



Thermally excited particles incident on the surface 
will be reflected unless they have an energy of at least 
2{VmA — e°)/(jh + 1). At low enough temperatures, 
not many particles will have enough energy to penetrate. 
This leads to a significant decrease of heat conductiv- 
ity through the interface, a phenomenon well-known for 
superconductors [isl . [ioj . We will now argue that the 
same blocking of energy transport occurs for imbalanced 
fermion gases, and discuss the influence of the interaction 
regime on this phenomenon. 

Heat transport across the interface can be studied us- 
ing the transmission coefficients calculated in last section. 
Using the relation: 



K 10° 

N 10 
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< Deep BCS 1/M=-2 




□ BCS with l/k^a =-0.44 
O Unitarity with 1/l<pa=0 
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FIG. 4: The thermal conductivity n across the N-SF interface 
divided by the normal-phase conductivity k,^ against k^T j ^ 
for the unitary, BCS and deep BCS regimes when m| = 
and V},i,s ~ 0. For k^T < O.IA, drops dramati- 

cally (notice the logarithmic scale). The curve represents the 
analytical result obtained using of the Andreev approxima- 
tion, Eq. (g3J. At unitarity, A = Omk^T^, in the BCS case 
A = 0.25kgTj, and in the deep BCS case A = 0.002fcgT^. 
For large values of the temperature, one finds that n/n^ = 1. 



one can write the particle-like heat flux through the in- 
terface and to the a branch as: 



7r2?l3(TO+l)2 



(21) 



where / is the Fermi-Dirac distribution and the integra- 
tion is performed over all regimes of the Xp~^a plane. 
For the heat flux carried by hole excitations, the in- 
tegration is only over regime VI and by use of dfc^ = dk^ , 
we can obtain the same expression as Eq. (j2ip but with 
instead of m| and with instead of S^, as given 
in (|31d[) . By the same approach, the contributions of the 
/3-branch are then found by performing the transforma- 
tions /ig — "T- ^ l/^ii a ^ P and t^i- 

In equilibrium, quasiparticles are incident on both 
sides of the interface such that the N-SF flux is balanced 
by an equal SF-N flux. Suppose, however, that the tem- 
perature on the N side is ST higher than on the SF side 
(and 5T/T < 1). A net heat flow SQ = + then 
results from the scattering to the a-channel {SQ" ) and to 
the /3-channel (SQi^). For small temperature differences, 
we may write JQ" — k" ST with k° the heat conductiv- 
ity caused by scattering to the a channel. Within the 
Kapitza approach, we have: 



d(W^+W;;;) 



(22) 



In the deep BCS regime, where the Andreev approxima- 
tion may be used, one may obtain an analytical expres- 



sion for K. Taking k'' = k*^ = k^ = /c,^ and assuming the 
temperature to satisfy k^T ^ A <C /i^ effects regime 
VI of Figs. [2] and [3] to be the sole important scattering 
regime; in other words, all reflections arc due to Andreev 
reflection. We obtain for the heat flux to the a channel 
(see Eq. (EH): 



-2/3(Avni-eO )/{fh+l) 



2Tr^ n^y/m + 1 
3A:^r(m + l)/4](fc^T)3/2 



For the a-channcl the heat conductivity of Eq. ([22|l in 
the BCS regime is then obtained: 



^/^^{■m^ + mi)e 



-2/3(A\/m-£0 )/(m+l) 



7r2fi3(m -1-1)3/2 

(AVm - e° )[A\/m - e° + 3k^T{fh + l)/2] 



(23) 



k^TAv m 



which amounts to Andreev's result when /i, = and 
r7i| — mi The heat flux and heat conductivity of the 
/5-channel is again found after the transformation fh — > 
1/to, t^i? fT-s ~^ ai^d a ^ f3. The energy carried 
by the /3 branch is a factor e~^'"'s/('™+i) lower than that 
of the a branch and due to coexistence /i ~ A and since 
A T, it can be neglected. 

Eq. ([23|) demonstrates that k decays exponentially fast 
with decreasing temperature. This can be understood as 
follows. For low enough temperature i.e. k^T <^ ^/mA, 
only incident particles (and holes) with energies slightly 
above the threshold ~ ^/mA can contribute to the 
heat conductivity; they have, however, a very low statis- 
tical weight: 



/ cx e 



-2/3(VmA-eO )/(m+l) 



< I. 



(24) 



It is this weight factor which also appears in expres- 
sion (I23p . The question now is: what will happen 
to this exponential decay if we tune the interactions 
from the BCS to the unitary regime? In order to an- 
swer this, we have numerically calculated the ratio of 
K = K° -f K^* to the conductivity in the N phase k^^ — 
kgT{fi-im'i + /Lt|m|)/7r2fi3, and show this in Fig. 3] as 
a function of k^T/A and we assumed m| = and 
^^T,i,S' — 0. We find that decreases drastically he- 

low k^T K 0.1 A and is independent of the interaction 
parameter k^a. Also in Fig. [5] we plot the heat conduc- 
tivity of Eq. (|23p: it is seen to give a good flt for all 
scattering regimes [gO]. The reason for this good agree- 
ment is that the threshold line = \/mA of the dia- 
gram of Fig. [2] is unaffected by changing the scattering 
regime. Furthermore, relation (|20p which describes the 
transmission coefflcient slightly above the threshold line 
guarantees that displays the same behavior in regimes 
V and VI. 
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FIG. 5: The thermal conductivity k across the N-SF interface 
divided by the normal-phase conductivity k„ against T /T^ 
at the interface in the unitary regime. The full line depict 
the Andreev approximation wherein the coexistence condi- 
tion is taken from the one-channel model. The dotted line is 
the Andreev approximation with the coexistence condition is 
taken from Monte Carlo (MC) simulation s 1231 and a super- 
fluid local-density approximation (SLDA) [33]. At the lowest 
experimental values, that is when T < O.OSTp, the conduc- 
tivity is already very low. 



As argued before, the main effect of varying l/{kpa) 
from large and negative (deep BCS) to zero (unitarity), 
is to shift the curves 1-4 in diagram of Fig. [^ Both the 
strong exponential decay in VmA/{k^T) and the square 
root dependence of on — VmA are unaffected by the 
variation of l/{kpa) and the introduction of the Hartree 
potentials, hence the very similar behavior of for 
all scattering regimes under study. On the other hand, 
what is different between the scattering regimes is the 
value of A; it is small in the sense that A <C ^bTV the 
BCS limit and is comparable with k^T^ at unitarity. 



V. EQUILIBRATION TIMESCALES 

The introduction of the Hartree potentials has two im- 
plications on the heat resistivity. First of all, the scat- 
tering diagram will change by horizontal shifts of lines 
1-4; as argued before, this does not affect the steep ex- 
ponential decay of resistivity. Secondly, the coexistence 
condition will change, that is, the relation between A 
and h will change. This, on the other hand may strongly 
affect the exponential decay as is seen from Eq. (|24|) . At 
unitarity, the relation between A and h is known from 
Monte Carlo (MC) [22| simulations and a superfluid lo- 
cal density approximation (SLDA) [33]. From these we 
adapt the values A = 0.504 /c^T^, = -0.516 fcr 
Ms = 0.936 /c^T^ [13 and finally = QA%{^J.^ + m)^. 
In Fig. [5l we compare the Andreev approximation with 
values of the one-channel model (full line) and of the MC 
and SLDA results (dashed line). We see that the one- 
channel model predicts the sharp decrease of k to be at 



higher temperatures than the MC and SLDA results. 

We find that the interface conductivity at the lowest 
experimental temperatures is very low in the sense that 
the resulting characteristic equilibration time is of order 
of seconds i.e. the trap lifetime. 

We explain now how we estimate the characteristic 
equilibration time r. Consider again a temperature bias 
5T across the interface. We are interested in the re- 
sulting total heat flux 5W , induced by the excess heat 
X]i=| J, ^Qi i'^ the N phase. By definition, the heat flux 
W is the amount of heat transferred per unit time per 
unit surface area. One can therefore estimate r by the 
relation 5W — Ylii^] i^Qili'^^) with A the interface 
area. Moreover, at position r in the degenerate N phase 
i =T,i, one has 6Q,{r) = t:"^ p^{T)T5TAT/[2^l,{r)] with 
Pi the particle density and i =1, | [5l| . Combining this 
with 5W = KkgST, one finds 

^^^!^ V f ^dr 

Assume now an isotropic trap of lO'' particles at polar- 
ization 0.5 with radial trapping frequency 458Hz such 
that the local density approximation is valid Q . For this 
configuration, we calculate that the characteristic time of 
relaxation r is of order seconds when T/T^ < 0.02. This 
means that according to our estimates, the relaxation 
time is of the order of the lifetime of the system at the 
lowest attainable experimental temperature T/T^ = 0.03 
such that in these experiments, the N and SF phases may 
not be equilibrated. 

At least two interesting experimental consequences of 
the low thermal conductivity for ultracold gases are dif- 
ferent to what is the case for a SC. First of all, in a 
superconductor, the heat conductivity has a lattice com- 
ponent which dominates the electron component at low 
temperatures; the lack of such a component in the sys- 
tem under study makes the decrease of the conductivity 
more significant [soj . Secondly, at very low temperatures, 
a N-SF temperature gradient will entail a pressure differ- 
ence 6P across the interface, and therefore the physical 
movement of the interface; this is impossible for normal- 
superconductor junctions. Taking 5P = (dP/dT)ST and 
an infinite thermal resistivity, one finds [5l| : 

'^^-^E^.o^x^f— )' (25) 
6 .f^^ Pi V M» / 

where P.„ is the T = pressure of species i =t: i- De- 
pending on the time scales involved, mechanical equilib- 
rium might be reached before thermal equilibrium. 

VI. DISCUSSION 

The step function model for the gap (see Eq. (Ilip ) is 
adequate for tackling the problem at hand, that is for 
estimating the heat conductivity across the interface. 
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First of all, this model should gives qualitatively cor- 
rect results as long as the quasiparticle wavelengths are 
of the same order or less as the length over which the gap 
varies. Indeed, this is clearly valid in the BCS limit, when 
the healing length for the gap, further denoted as ^, is of 
order = ti^kp /m^A and A ^ H'^k'^/m_^. Moreover, 
near unitarity, the step-function model is also justified 
since the length of variation of the gap is smaller than 
the quasiparticle wavelength [20|. Therefore, the step 
function of Eq. pT|) serves as a good approximation for 
all regimes under consideration. 

Second, the step-model for the gap is justified when- 
ever the interface width must be small as compared to 
the variation length of the trapping potential. This is 
also realized in experiments. 

We argue now that, even if we would take into ac- 
count the exact interface profile, our results would hardly 
be affected. By solving the Bogoliubov-de-Gennes equa- 
tions in a self-consistent manner, the interface profile for 
the gap could be determined [H^, H^- At zero tem- 
perature, the concomitant interface width is of order 

— h^kp/m^lS.. As was shown in Ref. [sj, such smooth 
interface gives rise to a zero transmission coefficients for 
particles with energy below the gap, and above the gap, 
smaller transmission coefficients than predicted by our 
step-function model. This simply means that particles 
incoming to the interface with energy below the gap can- 
not cross over to the superfluid phase and therefore they 
cannot carry energy across the interface. Hence our con- 
clusions of a blocking of energy transport across the in- 
terface remains valid, independent of the interface width. 
Note also that, due to the decreased transmission coeffi- 
cient above the gap, the heat conductivity will be even 
smaller than predicted here. 

The same conclusion can be drawn at finite tempera- 
tures: particles with energy below the gap cannot pene- 
trate the superfluid and the conclusions of Ref. [5J con- 
cerning the transmission coefficients will still be valid. In- 
deed, since it is known that the interface width increases 
with temperature as it varies between £(T = 0) = and 
C(T) (X 4(1 - T/T near T = T [21], i is certainly 

of the order of ^„ at the experimentally relevant temper- 
atures, that is below O.IT.. This can be seen in Refs. [stI 
and [5^ where the self-consistent calculations were per- 
formed for various temperatures and even for the strong 
coupling regime [57| . 

In our model, we do not take into account the tem- 
perature dependence of the gap. However, this is well 
justified as the gap A(T) depends very weakly on tem- 
perature at low T. In the BCS limit it is well-known that 
A(T) goes as A(T) - A(T = 0) cx -VTexp(-A(0)/T) 
(see for example in Ref. ISa ). Moreover, this weak varia- 
tion is also present in the strong coupling limit [11] and 
can be understood within the single-channel model from 
the appearance of the factor exp(— A/T) in the finite- 
temperature gap equation. 



VII. CONCLUSION 

We have studied the effects of a N-SF interface in ex- 
periments on trapped imbalanced fermion gases. We clar- 
ify the nature of the possible processes that may occur 
when a particle is incident onto the interface. Such par- 
ticles will be thermally excited at any finite temperature, 
and act as carriers of energy, thus eventually causing 
thermal equilibration between the N and SF phases. 

We find that reflection (both Andreev and normal, or 
specular) of low energy particles off the interface in ex- 
periments on trapped imbalanced fermion gases causes 
a suppression of energy transport from one phase to the 
other. At low enough temperatures (comparable to the 
currently accessible temperatures in experiments), this 
suppression grows exponentially with decreasing temper- 
ature. This, in turn, delays thermal equilibration of the 
system, which may result in a temperature difference be- 
tween the N and the SF. Our estimates of the timescales 
for this equilibration are in the seconds, which indicates 
that this effect is experimentally relevant. 

The incorporation of this temperature difference may 
allow current models to finally resolve the difficulties in 
explaining the experiments. 
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Appendix 

Assume particles of species t incident on the interface 
with momentum k^ and energy . Define first ^ = 
X±/A and: 

Jo = - k^){kp - kj ~ ,,ik- + k,){k- + fcj, 

J, = ^^{k- -k^){k-+k,)- ,,{k- -kj{k-+k^), 
J, = ,^{-k- - k^){k- +k,)- ,,{-k- - k,){k- + fcj. 

If the coordinate (Xp,£„) with Xp = fi^fc^/(2TO|) is posi- 
tioned in regime VI of Fig. [21 the matching of the wave 
functions and their derivatives at 2; = leads to the co- 
efficients: 



k,+ 
k,+ 

V 5k^,i J 



Jo 



/ ,^{k- + k^){k,-kp) \ 
' -.^(fc'.-f fcj(fc^-fcp) ^ 

—2mk (fc'' + fcp) 



(27) 
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For regime V (see Fig. [2]), we get: 



The transmission coefficients for the different regimes are: 



4:; 



J. 



For regime IV (see Fig. (2), we get 



/ ,Jk- + k^){-kp ~ kj \ 
-^^{k--k,){k^~kp) ^ 
-~2mk^{k'^ + fcf ) 

V -2,^k^ikp + kj J 



(28) 



J. 



' -^,X-k--K){K-kp) \ 

~2mk (— fc*" + fcp) 
\ ~2^Xik-+K) ) 



(29) 



For holes of species | which are incident on the interface, 
and which necessarily must have energy and momentum 
in regime VI, the scattering amplitudes are: 



Ap-" 

k.- 

Ap." 

k,+ 

V ^k^,i ) 



Bh.n 
k, — 



/ ^,(fc"-fcj(fc^ + fcj \ 
' -,Jk--k^)ikp + kJ ^ 

-2k^{k>^ + kp) 
-2k^ {k>^ - k ) 
V -2k,{kp+kl) J 



(30) 



For regime VI: = [4fcpfcp(fc'' + fc,J^(';f - m) (31a) 

+Ak-k^{k,~k^f{m~e^)\ ui. 

For regime V: Si = \\k^k^\k'^ - ^hl^l^^f - ™) (31b) 
+4fc'-fcjfc, +fc-|2(TO-<;f)]/|j;|2, 

^k^k^ik'- - k^f{m - (31d) 

+Ak-k,{k^ + k^f{e^-m)\ ui. 



For regime IV: 
For holes: 5'' 



Here we also used c — m/<, and < m < 
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